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Using measured events from the fluorescence detector of the Pierre Auger Obser- 
vatory, an unbiased distribution of the atmospheric slant depths where showers 
reach their maxima has been obtained. Analyzing the tail of this distribution the 
proton-air cross-section for particle production at center-of-mass energies per nu- 
cleon of 57 TeV is determined to be [505±22(stat)+28— 36(syst)] mb. Systematic 
uncertainties in the analysis arise from the limited knowledge of the primary mass 
composition, the need to use shower simulations and the selection of events. For 
the purpose of making comparisons with accelerator data we also calculate the 
inelastic and total proton-proton cross-sections using an extended Glauber model. 

1 The Observatory 

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic 
rays. The surface detector (SD), which detects secondary particles at ground level, 
consists of 1660 water-Cherenkov detectors spread over 3000 km 2 on a triangular 
grid with 1500 m spacing. The array is overlooked by 27 fluorescence telescopes (the 
fluorescence detector, FD) which measure the longitudinal development of extensive 
air showers (EAS) initiated by the ultra-high energy particles in the atmosphere 
above the array by detecting the fluorescence and Cherenkov light produced along 
the shower trajectory as the charged particles cross the atmosphere 1 . The FD 
can measure events with energies from approximately 10 17 ' 5 eV, while full trigger 
efficiency in hybrid mode (a fluorescence event in coincidence with at least one 
tank) is achieved at energies greater than 10 18 eV 0. 

2 Observables 

The number of particles in the EAS (or shower size) as a function of the atmospheric 
slant depth (the amount of atmosphere traversed from its upper edge in g/cm 2 ) 
is called the shower longitudinal profile. Most of the EAS energy is dissipated 
through the electromagnetic component. Therefore, the shower size increases until 
the average energy of the in the EAS is about the critical energy of 81 McV, 
when the rate of energy loss due to collisions and ionization begins to exceed that 
due to radiation. Photons have a similar attenuation length due to pair production. 
The slant depth at which the longitudinal profile of a shower reaches its maximum 
is called X max and is one of the primary observables of the fluorescence detector. 
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The distribution of observed values of A max is widely used to estimate the mass 
composition of the primary beam and it also carries information on the cross-section 
of the first interaction in the atmosphere |3I4) . 

The differences in A max between showers of the same primary energy can be 
due to fluctuations in interactions as well as to different masses of the primary 
particles. For purely proton primaries, the A max distribution is a convolution of 
the fluctuations in the shower development from the point of the first interaction to 
the shower maximum (which is dependent on details of the hadronic interactions) 
and the exponential distribution of the depth of the first interaction with mean free 
path Ap_ a i r = (ra a i r )/<7p_ a i r , with the mean target mass of air (m a i r ) ~ 14.5m p w 
24.3 b g/cm 2 and the proton-air cross section cr p _ a ir [4]. In a simple model, heavier 
nuclei can be considered as the superposition of many nucleons. For an iron nucleus 
with energy E, each nucleon would have an energy £7/56, and the superposition of 
56 lower energy subshowers would reach A max earlier (higher in the atmosphere). 

Because the shape of the X max distribution (at least in the range between 10 18 
and lO 18 5 eV which we consider) suggests that there is a significant fraction of 
protons in the primary beam (for details on the determination of the primary 
composition see [5]), we can obtain a proton-enriched sample by using only the 
deepest observed events. The remaining dependence on both the mass composition 
and possible photon contamination is evaluated using simulations and included in 
the systematic uncertainties. The slant depth range used in the analysis is defined 
by the slant depths of 20 % of the most penetrating events, which is a definition 
that can be easily applied also on simulated events. As the primary observable, we 
use the exponential slope A of the tail of the X max distribution [6] . 
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3 Event selection 

We use events collected by the fluorescence detector between 1 December 2004 and 
20 September 2010 that have a signal in at least one of the SD stations measured 
in coincidence with the FD. The geometry for these events is determined with an 
angular uncertainty of 0.6°. Further we reject events 

• without reliable measurement of the aerosol optical depth [7], 

• with excessive cloud coverage (more than 25 % of the sky), 

• when % 2 /Ndf is greater than 2.5 when the profile is fitted to a suitable function 
as this could indicate the presence of residual clouds, 

• with the total uncertainty of AT max greater than 40 g/cm 2 

• when the angle between the shower and the telescope is smaller than 20° due to 
the difficulties of reconstructing their geometry and high fraction of Cherenkov 
light in such showers, 

• when the reconstructed X max lies outside the field of view of the fluorescence 
telescopes. 

In total, applying these selection criteria yields 11,628 high-quality events [6]. 

4 Fiducial volume selection 

The fluorescence telescopes have a limited field of view in elevation ranging from 
about 2° to 30°. The geometrical acceptance as well as the limitations related 
to atmospheric light transmission introduce a bias into the distribution of X max . 
The bias is even enhanced by demanding that X max be within the observed profile. 
For example, many showers landing close to the FD will have their X max above 
the field of view. Some vertical showers, on the other hand, may have their X max 
below the ground. In both cases such events will be rejected and the observed X meLX 
distribution will be biased. 

To avoid such a bias, we apply the fiducial volume selection. First, we determine 
from the data the range of values where most (99.8 %) of the events lie (550-1004 
g/cm 2 ). Then, we select only showers with such geometries that each of the showers 
would be observable for any X max in the whole range of slant depths (irregardless 
of the actual X max of the shower). From this unbiased distribution of 1635 events 
we find that the 20 % most deeply penetrating events fall within the range 768 to 
1004 g/cm 2 . Repeating the same procedure for this range, we obtain a set of 3082 
events (shown in Fig. [T]), of which 783 events contribute to the estimated value of 
A, yielding 

A = 55.8 ± 2.3(stat) ± 1.6(syst) g/cm 2 , 

where the systematic error is estimated as the RMS of the distribution of A cal- 
culated using different selection procedures. The average energy of these events 
is io 18 ' 24±0 ' 005< ' stat ' ) eV, corresponding to a center-of-mass energy of 57 ± 0.3(stat) 
TeV in proton-proton collisions. 



4 



Jan Ebr for The Pierre Auger Collaboration 




800 



Figure 2. Relation between the slope A in the simulations and the proton-air particle production 
cross-section [8]. 
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Figure 3. Resulting proton-air particle production cross-section compared to other measurements 
and model predictions. The inner error bars are statistical, while the outer ones include also the 
systematic uncertainties for a helium fraction of 25 % and 10 mb for the systematic uncertainty 
connected with the fraction of photons. 



5 Conversion to cross-section 



We have used the four high-energy hadronic interaction models: QGSJET01 [5], 
QGSJETII.3 QTj], SIBYLL 2.1 [H , and EPOS1.99 [H] to perform Monte Carlo 
simulations, where we smoothly modify all hadronic cross sections by an energy- 
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Table 1. Summary of the systematic uncertainties 



Description Impact on 



A systematics 


±15 mb 


Hadronic interaction models 


-8, +19 mb 


Energy scale 


±7 mb 


Conversion of A to cTp™ air 


±7 mb 


Photons, < 0.5 % 


< +10 mb 


Helium, 10 % 


-12 mb 


Helium, 25 % 


-30 mb 


Helium, 50 % 


-80 mb 



Total (25% helium) -36 mb, +28 mb 



dependent factor /, given by 

In (£/10 15 eV) 



/ (E, / 19 ) - 1 + (/i 9 - 1) 



ln(10 19 eV/10 15 eV) 

(/ = 1 for E < 10 15 eV) where E denotes the shower energy and /19 is the factor 
by which the cross section is rescaled at 10 19 eV. For each hadronic interaction 
model, we find a value of f%g that reproduces the measured value of A and use 
it to evaluate the cross-section at E = 10 18 24 eV. The proton-air cross sections 
for particle production obtained by this procedure deviate from the original model 
predictions by less than 5 % in all cases, except for SIBYLL, for which it is 12 % 
smaller than the original prediction. 

Additional sources of systematic uncertainties are (see Table 1): 

• the systematic uncertainty of 22 % in the absolute value of the energy scale of 
the primary particles [13j . 

• the dependence of the procedure for retrieving the cross-section from A on 
additional parameters such as the energy distribution, energy and X max reso- 
lution, 

• the presence of photons (observational limits on the fraction of photons are 
< 0.5 %) and the unknown fraction of helium in the primary beam. 

The resulting value for the proton-air particle production cross-section is 

crP™ a d h . = 505 ± 22(stat)i 28 (syst) mb 
at a center-of-mass energy of 57 ± 0.3(stat) ± 6(syst) TeV. 

6 Comparison with accelerator data 

We also calculate the inelastic and total proton-proton cross sections using the 
Glauber model extended by a two-channel implementation of inelastic intermedi- 
ate states to account for diffraction dissociation Q3]. This calculation is model- 
dependent since neither the parameters nor the physical processes involved are 
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Figure 4. Comparison of the estimated inelastic pp cross-section, accelerator data and some model 
predictions. The inner error bars are statistical, while the outer ones include also the systematic 
uncertainties. For references to the data, see \E\. 



known accurately at the cosmic-ray energies. In particular, this applies to the elas- 
tic slope parameter B, the correlation of B with the cross section, and the cross 
section for diffractive dissociation. The result for the inelastic proton-proton cross 
section is 

(T^p = 92 ± 7(stat)±? 1 (syst) ± 7(Glauber) mb 

and the total proton-proton cross section is 

CT *°4 = 133 ± 13(stat)±^(syst) ± 16(Glauber) mb 

These results agree with extrapolation from data at LHC energies to 57 TeV for 
a limited set of models. We emphasize that the total theoretical uncertainty of 
converting the proton-air to a proton-proton cross-section may be larger than esti- 
mated here, because there are further possible extensions to the Glauber model to 
account for a variety of other effects, as discussed in [6]. 
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